Evolutionary Biology

The general random walk model (GRW) of Hunt (2006) is used to infer directional evolution in mean trait values from sparse fossil data by modeling phenotypic change as the accumulated result of small steps with mean step sizes and step variances. Using simulations and real data cases, Ergon (2026) showed that the step variances can be estimated reasonably well only when the mean trait values have small measurement errors, while for fossil data with realistic measurement errors they appear to be extremely difficult to find, and they are often found to be negative. In the simulations Ergon (2026) assumed that the true phenotypic mean values were known. Here, I essentially repeat these simulations under the assumption that only mean trait values with large measurement errors are known, and based on weighted mean squared error (WMSE) comparisons the conclusion is that weighted least squares (WLS) is a better method than GRW. A second conclusion is that WLS is a better method also in the possibly rare cases with large measurement errors where the GRW parameters are estimated well. The GRW method is simply not flexible enough to handle such cases. A third conclusion is that Akaike Information Criterion (AIC) results for GRW models with large measurement errors relative to the step variance may be overly optimistic.

The diversity of body shapes is one of the most prominent features of phenotypic variation in mammals. Yet, mammalian body shapes are poorly quantified and the underlying components contributing to its diversity as well as its relationship to other components of the skeleton are rarely tested. Here, we use lagomorphs (hares, rabbits and pikas) as a model system to (1) investigate which components of the skeleton contributed the most to body shape diversity, (2) examine the relationships between body shape and relative limb lengths, and (3) test how body size, ecotype, burrowing behavior, and locomotor mode influenced variation in lagomorph body shape and appendicular morphology. We quantified the body shape and functional proxies of the appendicular skeleton in 40 lagomorph species from osteological specimens held at museum collections. Using phylogenetic comparative methods, we found the relative length of the ribs and elongation or shortening of the thoracic and lumbar regions contributed the most to body shape evolution across lagomorphs. Second, we found that only leporids (hares and rabbits) exhibited a significant relationship between limb length and body shape, where more elongate species exhibit relatively shorter forelimbs and hindlimbs. Lastly, we found that models incorporating body size were the best predictors of lagomorph body shape and the majority of the appendicular traits, whereas models incorporating burrowing behavior and locomotor mode were largely poor fits. Broadly, these results indicate that larger lagomorphs tend to exhibit more robust body shapes with longer, more gracile forelimbs, whereas smaller lagomorphs tend to exhibit more elongate body shapes with shorter, more robust forelimbs. Overall, this work contributes to the growing understanding of mammalian body shape evolution and demonstrates the importance of not omitting body size in ecomorphological analyses.

Morphological measurements underpin a wide range of ecological and evolutionary research, yet the manual landmarking workflows on which most morphometric studies depend remain a persistent bottleneck that limits both the pace and scale of biological research. Machine learning offers compelling solutions, but most automated landmarking tools require substantial computational expertise, creating a gap between technical capability and practical adoption by biologists. Here, we present LizardMorph, an integrated machine learning pipeline and web-based interface for semi-automated anatomical landmark detection on biological images. LizardMorph couples a fine-tuned ML-Morph shape predictor with an accessible, browser-based interface that enables researchers to upload images, review automated landmark predictions, interactively correct outliers through point-and-click editing, and export results in standard morphometric formats--all without programming expertise or local software installation. Using dorsal X-ray radiographs of Anolis lizards with 34 anatomical landmarks as a proof-of-concept, we show that the ML-Morph model achieves high predictive accuracy, with landmarks on well-defined skeletal structures predicted with 100% accuracy within a 1 mm tolerance threshold. A controlled user study comparing LizardMorph against traditional manual landmarking (TpsDig2) demonstrated significant efficiency gains: experienced annotators completed LizardMorph landmark verification 37.5% faster than manual annotation. Extrapolated to batch processing 1,000 lizards, LizardMorph saves experienced researchers approximately 6.5 hours of manual processing time. Critically, LizardMorph implements a human-in-the-loop design in which automated predictions serve as editable starting points, preserving researcher oversight and enabling correction of the occasional large-error outliers that would be unacceptable in fully automated workflows. LizardMorph is freely available as an open-source tool and provides a replicable framework for developing ML-assisted annotation tools that can democratize access to high-quality morphometric analysis across diverse biological research communities.

Exaggerated insect traits often show positive allometry, yet nonlinear scaling can reflect either continuous curvature or discrete morphs. Distinguishing between these alternatives is important because they imply different developmental and evolutionary scenarios. Using cross-sectional data from 1,000 adult Oryctes rhinoceros, we analyzed the static allometry of nine traits with pronotum width as the primary body-size proxy and body length for sensitivity analyses. Cross-validated comparisons among linear, continuous nonlinear, and two-component mixture models showed that continuous nonlinear models improved predictive performance over linear models for horn length and pronotal depression width in both sexes. By contrast, mixture regression did not outperform the best continuous model in either sex, providing no positive support for discrete within-sex dimorphism under the tested model set. After accounting for body size, horn length remained positively associated with pronotal depression width, indicating size-independent covariation. These associations were retained when body length was used instead of pronotum width, supporting robustness to body-size proxy choice. Together, these results support continuous nonlinear adult scaling of the horn and pronotal depression in O. rhinoceros and indicate covariation not attributable solely to body size under the tested model set.

Given the ever-increasing availability of highly detailed phenotypes, modelling trait evolution in a multivariate framework is becoming a challenging task. Current phylogenetic comparative methods often struggle with high-dimensional datasets because they suffer computational limitations and interpretability. Here, we propose a maximum likelihood-based approach called Probabilistic and Phylogenetic Principal Components Analysis (P3CA) to circumvent current limitations. This approach is based on a continuous latent variable model, whereby observed traits are explained by a smaller number of unobserved variables that evolve according to a given evolutionary model. We implement the approach under Pagels lambda model using an Expectation-Maximisation algorithm that makes it computationally efficient and allows missing values. Using simulations, we demonstrate that evolutionary parameters are accurately estimated, regardless of phylogenetic signal, the number of traits or the proportion of missing values. The reconstruction of the reduced space is more accurate than the one obtained using other dimensionality reduction approaches, such as phylogenetic and conventional PCA. Likewise, the estimated values for missing data are more accurate than those obtained using current phylogenetic data imputation approaches. We illustrate the approach on a 3D geometric morphometric dataset describing Crocodyliformes skull shapes and containing around 4% of missing data. Our P3CA method unlocks the possibility to analyse and more easily interpret the large-scale multivariate datasets generated in recent decades within a phylogenetic comparative framework.

Human dental morphology is diverse and varies both within and between populations worldwide. Common variants include different numbers of cusps and roots, as well as different configurations in the fissures, ridges, and grooves on tooth crowns. Because teeth preserve well in taphonomic contexts and retain strong genetic signatures in their morphology, researchers across disciplines use dental form for research ranging from population affinity identification in forensic cases to the reconstruction of population history in archaeological and paleontological studies, and to explore the genetic underpinnings of dental development. However, these analyses are limited, as no publicly available comprehensive database of human dental morphological variation currently exists. Data are typically shared only within professional networks, excluding scientists outside those circles. Information is dispersed across repositories and the supplementary materials that accompany peer-reviewed publications, creating a fragmented and difficult-to-navigate data landscape. Although numberous dental anthropologists devoted their careers to collecting extensive datasets from thousands of individuals worldwide, their data have not yet been published in raw form nor made compatible. Here, we introduce the Dental Morphological Database (DeMoDa), an open access repository comprising 246 dental traits for 17,308 individuals across 32 major geographic regions worldwide spanning the past several thousand years. These data are from the legacy datasets of Turner, Hanihara, Scott, and Irish. Because these researchers used different scoring systems, we developed a scoring system with two steps that enables the three most-common scoring systems to be combined. We provide two versions of the data that correspond to these two steps: one with more nuanced trait scores at the cost of a smaller sample size, and another with less nuanced trait scores but broader sample coverage. We discuss the implications of publishing these legacy data and outline our decision-making process that guided their release in accordance with both the FAIR and CARE principles of open science.

Mate preferences are often influenced by the magnitude of sexual signals, which are presumed to indicate underlying aspects of signaller quality. Although the perception of these signals depends on sensory processes, the role of perceptual adaptations and constraints in mate assessment is frequently overlooked. Many sensory systems follow Webers law of proportional processing, where discrimination between signals is based upon their proportional, or relative, difference rather than their absolute difference. Because preference strength varies with relative trait magnitude, Webers law could strongly influence sexual selection, changing the coevolution of traits and preferences. Here, we explore the consequences of Webers law for sexual selection using individual-based models, applying Scalar Utility Theory to mate choice. We investigate the coevolution of male ornaments and female preferences under both Fisherian and good genes scenarios, as well as scrutinizing the sexual selection of multiple ornaments and preferences. Including Webers law in these models either reduced ornament exaggeration, or promoted exaggeration and diversification of ornaments and preferences, depending on the costs of choice and how rapidly female survival decreases when preferences evolve away from the naturally selected optimum. These results highlight the importance of perception and cognitive processing in shaping sexual selection and its evolutionary impacts.

Understanding how anatomical structures evolve requires disentangling the roles of integration and modularity in shaping morphological variation. The vertebral column, a serially repeated and regionally differentiated structure, provides a powerful system for investigating these processes. Here, we examine how vertebral morphology evolves in relation to whole-body elongation across the adaptive radiation of Lake Malawi cichlid fishes. We tested for evolutionary integration between the precaudal and caudal domains, as well as assessed the contributions of vertebral count, centrum shape, and intervertebral spacing on body elongation. We find strong evolutionary integration between precaudal and caudal vertebral shape, with both vertebral shapes varying along shared axes of multivariate shape change. Despite this, precaudal and caudal vertebral counts evolve independently, indicating a decoupling between the evolution of identity and morphology. Whole-body elongation is significantly associated with coordinated changes in vertebral and rib morphology, including proportional increases in centrum size, posterior displacement of neural and haemal spines, and increased rib curvature. In contrast, centrum elongation and intervertebral spacing do not independently explain body elongation beyond vertebral counts. These results demonstrate that body elongation in cichlids necessitates integrated, multivariate changes in axial morphology. Our findings highlight the importance of morphological integration in facilitating coordinated evolutionary responses in anatomical systems.

Anseriformes (waterfowl) and Galliformes (landfowl) are among the worlds most recognisable groups of birds, together comprising the clade Galloanserae. Despite their close evolutionary relationship, the skulls of adult anseriforms and galliforms exhibit strikingly distinct morphologies, the developmental basis and evolutionary history of which is poorly understood. To illuminate the developmental and evolutionary underpinnings of cranial disparity between and within these major extant bird clades, we quantitatively investigated ontogenetic changes in cranial morphology across galloanseran phylogenetic diversity, focusing on the previously unexplored post-hatching interval during which adult morphology takes shape. Our results reveal the combined effects of multiple heterochronic shifts early in galloanseran evolutionary history including anseriform hypermorphosis, along with influential non-heterochronic changes leading to substantially more disparate ontogenetic trajectories--and greater cranial variability--in anseriforms than galliforms. Key galloanseran fossils help clarify the polarity of evolutionary shifts in cranial development through galloanseran phylogenetic history and demonstrate that extant galliform cranial morphology is more constrained and retains a more plesiomorphic morphology than that of anseriforms. Our work helps illuminate the developmental basis of the iconic differences in cranial form between waterfowl and landfowl and illustrates the importance of broad phylogenetic and ontogenetic sampling for clarifying patterns of post-hatching developmental divergence among major vertebrate clades.

Birds and mammals are shrinking and shapeshifting as global temperatures rise. Ecogeographic rules predict that such changes should ease heat stress by increasing surface-area-to-volume ratios, and thus, the capacity for heat exchange. This has led to the hypothesis that body size reductions are driven by thermoregulatory selection or adaptive plasticity, although recent syntheses point to more complex, multifactorial causes. Crucially, recent theoretical models predict that thermoregulatory benefits of smaller body size only emerge at extreme deviations from average phenotypes. Here, we exploit agricultural selection in Japanese quail to directly test this hypothesis, using three breeds spanning extreme differences in body mass, surface area, and relative appendage lengths. Evaporative cooling capacity and the scope for evaporative water loss broadly followed allometric predictions when contrasting small and larger breeds. As expected, this allowed the smallest breed to tolerate higher air temperatures. However, differences in heat tolerance limits between breeds were consistently much smaller than predicted. Additionally, the breadth of thermoneutral zones overlapped in full, and upper critical temperatures were remarkably similar, between breeds. Together, these results show that heat tolerance is only weakly linked to surface-area-to-volume relationships and cannot be explained by size alone. Thus, although smaller bodies may modestly enhance heat dissipation when size variation in a population is substantial, our findings suggest that recent body size reductions and morphological shifts are unlikely to be driven primarily by thermoregulatory benefits.

Convergent evolution is often interpreted as evidence that similar ecological challenges select for optimal functional solutions, yet the biomechanical mechanisms underlying such convergence remain poorly resolved. Parrotfishes independently evolved fused, beak-like dentition multiple times, a trait associated with feeding on hard substrates such as coral, but its advantage has not been quantified. Here, we combine finite element analysis with shockwave modeling to evaluate the performance of the beak during feeding. Across five species spanning independent evolutionary origins of beaked and non-beaked dentitions, we find that jaw opening is governed primarily by geometric scaling, with minimal influence from dentition morphology. In contrast, under stabilized biting, fused dentition reduces the rate of stress accumulation relative to discrete teeth, indicating a functional advantage under constrained loading. Shockwave modeling further shows that tooth geometry and stacking regulate stress propagation: smooth profiles reduce stress concentration, while stacked architectures localize stresses within the tooth material and limit transmission into the surrounding bone. These results suggest that convergent beak evolution in parrotfish reflects repeated optimization for stress management under constrained and impact loading, and reveals a general principle by which biological structures control internal force transmission during high-force feeding rather than maximize strength alone.

O_LIAutomated landmarking transfers anatomical landmarks from a reference specimen onto many targets, greatly increasing analytical throughput. However, this procedure needs to be bootstrap using an initial sample. An arbitrary or atypical choice imprints a reference-of-origin bias that propagates through the pseudo-landmarks, the resulting morphospace, and the downstream template selection, a risk that is difficult to avoid for large datasets whose variation is not yet understood. C_LIO_LIWe replace the fixed reference with an iterative consensus atlas, warped over a few iterations toward the Procrustes mean shape of all similarity-aligned specimens. We evaluated it on a 62-strain Mus musculus skull panel by running both the original fixed-reference pipeline and the new consensus pipeline 62 times each, using every specimen in turn as the bootstrap. We compared atlas convergence, inter-atlas similarity, morphospace reproducibility, reference-choice variance of pairwise Procrustes distances, downstream k-means selection stability, and leave-one-out out-of-sample fit, and tested generalisation on great-ape datasets of differing sampling balance. C_LIO_LIThe consensus atlas converged within a few iterations and was far less sensitive to the starting specimen than the fixed reference. It produced more reproducible morphospaces (mean RV 0.960 versus 0.944), reduced the reference-of-origin variance of pairwise distances by a median of about 60%, drew downstream template selections from a smaller and more consistent pool of specimens, and fit held-out specimens more closely in all 62 strains. On the great-ape data the atlases agreed closely in surface geometry, but the downstream morphospace became reference-dependent when the sample was taxonomically imbalanced, and a smaller balanced subset outperformed the larger imbalanced one. C_LIO_LIIterative consensus atlas building removes a persistent bias from automated landmarking and yields reference-invariant, reproducible results, with sampling balance mattering more than absolute sample size. Because the atlas stabilises quickly, it can be built from a small balanced subset while the remaining specimens are simply landmarked against it, a practical route to scaling reference-invariant landmarking. The method is implemented in ALPACA within SlicerMorph, with a mock library enabling headless use on HPC. C_LI

Cities are characterized by elevated temperatures, increased pollution, and high-density human populations which often are accompanied by changes in available resources, like food. These shifts have the potential to drive phenotypic divergence in urban wildlife. Functional morphological traits, like body size, can mediate interactions between wildlife and habitat and are closely tied to life history and fitness. While examples of functional morphological variation associated with urbanization are increasing, variation in such traits as a response to urbanization remains unexplored for most taxa. Here, we investigated morphological divergence between urban and rural populations of house mice (Mus musculus domesticus). House mice are globally distributed in diverse habitats and are a model system with a wealth of phenotypic data, making them useful for the study of the impacts of urbanization on morphology. Using a paired replicate design, we sampled urban and rural populations in three distinct metropolitan regions in the eastern United States. We found that body size was smaller in urban populations. Using 3D geometric morphometrics, we also analyzed variation in cranial shape across habitats. Differences in cranial shape were largely allometric, that is, driven by differences in body size. However, we also uncovered evidence of cranial shape variation between habitats not explained by size. In contrast, we did not find evidence for habitat-driven differences in cranial capacity independent of size. Overall, our results suggest a key role for body size in mediating morphological responses to urbanization and highlight the potential of house mice as a globally-distributed model for urbanization.

Within species, male testosterone is often linked to mating competition and paternal care, suggesting that sex differences in endogenous testosterone values across mammals may covary with broader reproductive strategies. Using a structured literature search, we compiled 63 studies, spanning 31 non-human species and 9 human populations, reporting endogenous, non-experimentally manipulated testosterone values for both adult males and females within the same population and context. From these studies, we calculated male-to-female testosterone ratios, and analysed these data using Bayesian phylogenetic multilevel models. We tested whether testosterone dimorphism was associated with paternal care and sexual size dimorphism while accounting for sampling matrix, assay method, breeding context, and wild versus captive setting. Across non-human mammals, neither paternal care nor sexual size dimorphism (indexing competition) showed a clear association with testosterone ratios, and the same pattern emerged in the primate-only subset. By contrast, sampling matrix was consistently associated with testosterone dimorphism across all analyses, with lower male-to-female ratios in non-blood than in blood-based measures. In primates, testosterone ratios were also lower in captive than in wild populations, although this pattern was not clearly supported in the broader non-human dataset. In the human-only analysis, testosterone ratios did not clearly differ between industrialized and small-scale societies, whereas the matrix effect remained evident. Overall, our results suggest that sampling matrix is a major source of variation even for ratio-based measures, highlighting the need for caution when inferring between-species endocrine differences from studies using different substrates. More broadly, directly comparable, non-experimentally manipulated testosterone data for both sexes remain rare across mammals, limiting comparative inference.

Although it is well established that certain stages of development are molecularly more conserved than others, the reasons for this phenomenon remain largely unknown. We study molecular conservation in the development of an organ, the molar, by comparing the temporal profiles of expression in mice and hamsters. We find that the cause of conservation of expression and of coding sequences changes over molar development. Gene expression levels display a classical increase of divergence as development progresses. In terms of genes expressed, the composition of early and late stages is better conserved and enriched in pleiotropic genes, yet each stage mobilizes different sets of pleiotropic genes, cell division for bud growth and secretion for tooth mineralization. Moreover similar patterns of higher divergence of gene sets and of coding sequences at mid development, are caused by different biological phenomena, in that case heterochronies and blood colonisation respectively. In conclusion, the patterns of molecular conservation in developing molars are shaped by a combination of processes intrinsic to the teeth, and by negative and positive selection on functions which are mostly extrinsic to the teeth. This is likely translatable to explain molecular conservation patterns in many other biological systems. AUTHOR SUMMARYFor species to evolve different adaptations to different life styles, their anatomy has to evolve correspondingly. This in turn implies evolution of the embryonic development of anatomical structures. Notably, tooth shape can evolve rapidly as an adaptation to different diets. Mice and hamsters are closely related rodents who yet differ in the shape of their molars, and thus in their development. In this study, we investigated why the genes active in molar development are more or less similar between the two species from early tooth bud to fully formed embryo molar. We found that early and late molar development were slow evolving, while mid-development was evolving faster. But surprisingly, this was in part due not to tooth evolution, but to the involvement of genes which are active in other processes in the body. For example an influx of immune cells also brings fast evolving immune genes. This helps us understand better the complexity of causes of apparently simple evolutionary patterns.

Cephalo-pelvic disproportion in humans has traditionally been interpreted through the obstetrical dilemma framework, assuming a trade-off between bipedal locomotion and childbirth. However, cephalo-pelvic covariation and pelvic sexual dimorphism might be common adaptations to parturition among mammals. We use a controlled hybridization model in mice to test whether cephalo-pelvic covariation and pelvic sexual dimorphism are population-specific, genetically structured, and sensitive to hybridization. We analyzed skull-pelvis variation and covariation, as well as sexual dimorphism of pelvic morphology across four divergent wild-derived mouse strains and their hybrids. Hybridization induced consistent cranial and pelvic size enlargement. Females exhibited significant cephalo-pelvic shape covariation, characterized by an association between rounder, wider birth canals and larger neurocrania, consistent with functional integration under obstetric selection. Hybrids showed disrupted size covariation, increased pelvis shape variance, and reduced cephalo-pelvic integration. Pelvic sexual dimorphism was systematically reduced in hybrids. Cephalo-pelvic covariation and pelvic sexual dimorphism are not exclusive to bipedal or encephalized species. They likely reflect widespread selection on birth canal morphology in mammals and have a complex genetic basis sensitive to hybridization. These findings weaken a human-exclusive interpretation of the obstetrical dilemma and highlight genetic introgression as an understudied factor shaping cephalo-pelvic integration and disproportion risk in mammals, including humans.

Phylogenetic comparative methods are a critical tool in biology, providing the framework to test evolutionary hypotheses of phenotypic diversification. Accommodating intraspecific variation in these analyses is critical for accurate evolutionary inference, but current multivariate methods either assume traits evolve independently or that all taxa share the same intraspecific covariance structure. Violations of these assumptions can produce biased estimates of evolutionary parameters. Here, we introduce a hierarchical Bayesian framework for multivariate traits that jointly estimates taxon-specific intraspecific covariance structures alongside the underlying evolutionary process. This framework propagates uncertainty from sample size discrepancies and missing data, enabling the incorporation of highly variable morphological traits into phylogenetic analyses. Analyses of simulated data demonstrate that our framework achieves well-calibrated coverage (95%), whereas the standard practice of treating taxon means as known without error reduces coverage to 70%, confirming that ignoring intraspecific variation produces systematically biased evolutionary inference. We apply this framework to estimate the evolutionary rates and intraspecific distributions underlying perikymata distribution diversity in great apes, including modern humans and Neandertals. We show that, compared to other great apes, canine perikymata spacing in the genus Homo likely evolved under a substantially different regime than non-human apes, with cervical enamel evolving both rapidly and in a coordinated, modular fashion. We further find that Gorilla and Pongo show striking conservation in perikymata spacing relative to late Homo and Pan. These results and our method, which is applicable to other multivariate traits, provide the first phylogenetically rigorous characterization of an enamel growth trait across the great ape clade and establish intraspecific uncertainty propagation as a necessary component of multivariate phylogenetic analysis.

Accurate biological sex estimation is a key objective in archaeological and bioanthropological research but remains challenging when skeletal remains are fragmented, juvenile, or poorly preserved. Paleoproteomics approaches based on the detection of sex-specific amelogenin peptides (AMELX/AMELY) have emerged as a powerful alternative to osteological and genetic methods. However, current workflows often lack standardized criteria for peptide-level confidence assessment, potentially affecting the reproducibility and reliability of sex assignments. In this study, I evaluated the impact of peptide-level confidence filtering on paleoproteomics-based sex estimation through the reanalysis of 164 Homo sapiens individuals from 10 published datasets and 26 Bos taurus individuals from 3 datasets, spanning contexts from the Pleistocene to the present. To address methodological inconsistencies, I developed SexPeptID, an R/Shiny-based framework that integrates Posterior Error Probability (PEP) filtering, standardized peptide selection, and explicit uncertainty assessment. Application of SexPeptID revealed that peptide-level filtering substantially affects sex assignment outcomes: 17 previously classified males (10.4%) were reclassified as non-conclusive, while 5 individuals (3.1%) were identified as potentially female. Despite this sensitivity, AMELX/AMELY-based sex estimation remained robust overall, with stable signal ratios observed across archaeological periods. Variability in peptide intensities was primarily associated with dataset-specific factors rather than temporal differences, highlighting the influence of analytical workflows and preservation conditions. By incorporating confidence-based filtering and a non-conclusive classification category, SexPeptID improves the transparency, reproducibility, and reliability of palaeoproteomics sex estimation, providing a standardized framework for future archaeological and bioanthropological studies. HighlightsO_LISexPeptID provides a reproducible framework for amelogenin-based sex estimation. C_LIO_LIPeptide-level confidence filtering significantly affects paleoproteomics sex estimates. C_LIO_LI13.4% of published male assignments were revised after confidence filtering. C_LIO_LIAMELX/AMELY ratios show temporal stability from modern to Pleistocene samples. C_LIO_LIStandardized uncertainty assessment strengthens palaeoproteomics inference. C_LI

Sharks exhibit extraordinary morphological diversity across a wide range of ecological niches, yet large-scale, high-resolution digital datasets of their internal anatomy remain limited. Here we present an open-access 3D shark anatomical repository derived from published X-ray computed tomography (CT) data, featuring manually segmented and systematically annotated models of the chondrocranium, visceral arches, axial skeleton, musculature, and viscera in standard STL format. The dataset comprises 117 individuals, representing 72 species across 25 families and all nine extant shark orders, with 115 full-body reconstructions and two head-only models. This open-access dataset offers a comprehensive resource for comparative anatomy, biomechanical simulations, evolutionary developmental biology and biomimetics research of extant sharks.

Matching habitat choice provides a mechanism for individuals to maximise their expected fitness by selecting an environment that better fits their phenotype. Many animals choose their local environment by evaluating levels of perceived predation risk against possible resource gain. To test if predation risk is a major driver of habitat choice, we quantify scototaxis, or preference for dark versus light backgrounds, in juvenile guppies. As light backgrounds increase visibility to predators, this aspect of habitat choice captures variation in boldness in small fishes. By rearing and testing 586 fish descended from ten natural populations from Trinidad under common garden conditions, we first quantify (broad sense) heritable variation, i.e. evolutionary potential, within populations. Next, we test for evolutionary divergence among populations in mean preference, and if present, whether ancestral predation regime is a mediator of divergence. Finally, we ask whether families and/or populations differ in the amount of behavioural variation they contain. Habitat choice varied among families (12% of total variance), consistent with heritable variation (0.2). We also found mean preference varies among populations (11% of total variance explained). Evolutionary divergence among-populations is partly explained by ancestral predation regime, with populations from low-predation sites showing a stronger average preference for dark backgrounds than high-predation populations from the same river. Additionally, we find that within-population behavioural variation is greater in high-predation populations. We conclude that guppy populations contain heritable variation that could facilitate adaptive evolution if scototaxis is subject to natural selection. Furthermore, while genetic drift may also contribute to evolutionary divergence among-populations, observed patterns are qualitatively consistent with local adaption to predation regime. Our results suggests that high predation sites favour bolder habitat choice on average, but also that local predation regime shape the evolutionary dynamics of variation, perhaps by maintaining shy-bold variation among-individuals or by favouring individuals with less-predicable behaviour.